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Abstract 

The end stage of a gravitational collapse process can generically result in 
a black hole or a naked singularity. Here we undertake a comparative 
analysis of the nature of tidal forces in these backgrounds. The effect of 
such forces is generically exemplified by the Roche limit, which predicts 
the distance within which a celestial object disintegrates due to the tidal 
effects of a second more massive object. In this paper, using Fermi normal 
coordinates, we numerically compute the Roche limit for a class of non¬ 
rotating naked singularity backgrounds, and compare them with known 
results for Schwarzschild black holes. Our analysis indicates that there 
might be substantially large deviations in the magnitudes of tidal forces 
in naked singularity backgrounds, compared to the black hole cases. If 
observationally established, these can prove to be an effective indicator of 
the nature of the singularity at a galactic centre. 
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1 Introduction 


Einstein’s theory of General Relativity (GR) [T] is the most successful theory of 
gravity till date. An intriguing prediction of the theory is the existence of math¬ 
ematical singularities, which have been a major focus of research in gravity over 
the last century. Indeed, black holes, which are singular solutions of the Ein¬ 
stein’s equations predicted by GR, form a basic ingredient in our understanding 
of the universe as it is generally believed that the centre of supermassive galaxies 
are candidate black holes. There are also other singular solutions of GR called 
naked singularities that are not as well understood, but their existence cannot be 
ruled out. Possible distinctions between black holes and naked singularities have 
been of great interest in the recent past. To substantiate and carry forward such 
analysis is the purpose of this paper. 

Recall that singularities in GR arise in gravitational collapse processes, see, 
e.g |2]. Indeed, In the absence of pressure, a spherically symmetric distribution 
of matter will collapse under its own gravity. This scenario is often repeated in 
more realistic situations including effects of pressure. The end stage of a generic 
collapse process is singular, and results in a black hole or a naked singularity [3] , 
depending on initial conditions. While black holes are space-time singularities in 
GR that are covered by an event horizon, naked singularities do not have such 
a cover. The cosmic censorship conjecture originally proposed by Penrose more 
than four decades back states that nature does not allow naked singularities. 
However, a complete proof for this conjecture is lacking till date, and several 
research groups have directed their attention to the formation of naked singu¬ 
larities, their stability, and in particular distinguishing features between naked 
singularities and black holes. This last points assumes signihcance in predict¬ 
ing observational differences, if any, between black holes and naked singularities. 
Such analyses has been carried out for example in the context of accretion discs |1] 
and gravitational lensing 0 . 0 . In this work, we focus on another important 
aspect of gravity, namely tidal forces, and the purpose of this paper is to investi¬ 
gate the difference in the nature of these forces in naked singularity backgrounds, 
compared to black hole cases. The motivation for this line of approach is that 
such differences might be signihcant to be observable. In particular, tidal effects 
on neutron stars [7] in black hole backgrounds have been extensively studied in 
the literature and one of our aims is to quantify the differences between these vis a 
vis naked singularity backgrounds, which might be relevant in realistic situations. 

Tidal forces are manifestations of non-local gravitational interactions. Gon- 
sider, for example, a celestial object moving under the gravity of a more massive 
object. Non local gravity might cause this object to disintegrate. Viewed in the 
Newtonian context, this is easy to visualize. Gonsider for example a massive star 
of mass M and radius R that has a satellite made of incompressible matter, of 
mass Mg and radius rj,, both objects assumed to be spherical, with their centres 
at a distance d apart. An estimate of the tidal force can be obtained by comput¬ 
ing the gravitational forces due to the star at the (near) end of the satellite and 
its centre. Assuming that Tg/R -C 1 and equating this with the self gravitational 
force of the satellite predicts that the satellite will disintegrate if the dimension¬ 
less ratio {Mg/M){d/Rg)^ is less than a critical value, this value depending on the 
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Figure 1: Critical Mass For Circular geodesics in the Kerr geometry for different 
values of a, the rotation parameter. The singularity at r = 0 is naked for a > 1. 


spin angular velocity of the satellite. Then, this ratio translates into the fact that 
tidal disruptions occur if the radial separation between the star and the satellite 
is less than a critical value, known as the Roche limit. 

While tidal forces are easy to understand in the Newtonian context when the 
stellar mass is made of incompressible fluids, they may be substantially difficult 
to compute in the framework of GR, especially if one considers the fluid dynamics 
of the stellar object. The importance of this latter fact has been recognized of 
late, and by now a large body of literature is available, which considers the nature 
of the fluid star that is being tidally disrupted. Appropriate coordinates to deal 
with the problem in the framework of GR are the Fermi normal coordinates [8]. 
Here, one hrst sets up a locally flat system of coordinates along a geodesic, and 
then computes the tidal force taking into account the internal fluid dynamics of 
the star, and hence obtain the Roche limit. The actual computation procedure 
is non-trivial, and with little analytical handle available, one usually resorts to 
numerical analysis, along with a number of simplifying assumptions. This method 
was hrst elaborated upon in the work of Ishii et al [9], where tidal effects on objects 
in circular geodesics around a Kerr black hole were computed, upto fourth order 
in the Fermi normal coordinates. This built upon the work of Fishbone [lO], and 
to the best of our knowledge, was the hrst GR computation of tidal ehects taking 
into account the equilibrium huid dynamics of the stellar mass. The huid star 
was treated in a Newtonian approximation, wherein it was possible to superpose 
a tidal potential calculated in Fermi normal coordinates on the star’s Newtonian 
potential, and a numerical routine was used to extract various physical quantities 
related to the tidal force. 

The purpose of this paper is to generalize and extend the analysis of [9] to 
naked singularity backgrounds. A simple possibility for example is to consider 
tidal forces in the Kerr background considered in [9], and extend it to the regime 
where the singularity is naked. This might happen for example, when a black 
hole accretes angular momentum by capturing rapidly spinning stellar objects. 
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High spin black holes are known to exist in nature, and studies involving the spin 
of the Kerr black hole lesser than, but close to its mass, have appeared in un- 
Here we investigate the situation in which the spin of the black hole exceeds its 
mass. 

We will, for most of this paper, adopt geometrized units in which the Newton’s 
constant and the speed of light is set to unity. Then, the Kerr solution given in 
Boyer Lindquist coordinates reads 


ds^ = - 1 - 


2Mr\ 


P" 


■) 


dt^ + ^dr'^ + 
A 


(r^ + — Aa^sin^^ . 




sin^Odcj)^ 


n ,„n AMrasin^O , , , 

+ p'^dO'^ --- dtdcj) , 


( 1 ) 


with = r^ + a^cos^^, and A = r^ + a^ —2Mr. Here, M is the ADM mass and a, 
the angular momentum per unit mass of the source is a parameter that denotes 
the spin of the black hole. For a > M, the singularity at r = 0 is naked, and there 
is no event horizon. We will also set the ADM mass of the singularity, M = 1 
and the radius of the fluid star to 0.5, a choice made consistently in all examples 
considered in this paper. Fig.(IT]) (obtained by methods illustrated in the later 
sections where we have used the Fermi normal coordinates upto second order) 
shows the variation of the critical mass below which a star in a circular geodesic 
orbit is tidally disrupted, as a function of the rotation parameter (circular orbits 
in the Kerr naked singularity background and their stability has been studied 
extensively in [12] . We use the result of that paper that all circular orbits in the 
Kerr naked singularity background are stable). Here we have chosen the radius 
of the orbit of the star as r = 8 (in geometrized units), and used a polytropic 
equation of state for the star (as given later in Eq. fl25l) ). The solid red, dot dashed 
blue and dotted green curves correspond to different values of the polytropic index 
(0.5, 1, and 1.5 respectively). The Kerr singularity is naked for a > 1. We see 
from the hgure that there is no special behavior of the critical mass near the region 
a = 1, and that the tidal force decreases with increasing angular momentum of the 
source. This latter fact implies that for the Kerr naked singularity background, 
tidal effects are smaller than those in Kerr black hole backgrounds. That this is 
not generically the case is one of the main results of this paper. 

Two things need to be kept in mind at this stage. First of all, we are using a 
vacuum solution of GR, which is an good approximation in the above example, 
with r = 8, i.e far from the central singularity. More realistic situations arise when 
there is diffused matter (possibly dark matter) around the central singularity, and 
this will be considered later in this paper. Secondly, as can be seen from Fig.([T|), 
the change in the critical mass as a function of the rotation parameter is small. 
Hence, there is little hope to distinguish between a Kerr black hole and a Kerr 
naked singularity from the point of view of tidal disruptions. As we will see later 
in this paper, other naked singularity backgrounds offer a better scenario, and 
that one could in principle obtain effects of tidal disruption that are orders of 
magnitude apart from corresponding black hole situations. 

This paper is organized as follows. In the next section, we will hrst summarize 
the main ingredients used for our calculations and the various space-times con¬ 
sidered in this work. This involves three main steps. First we set up the Fermi 
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normal coordinates for the naked singularity backgrounds. Next, the hydrostatic 
equations for a fluid star at equilibrium are set up, and in the final step, a nu¬ 
merical routine is elaborated upon, for both circular and radial geodesic motions. 
The last two analyses discussed here essentially follow the work of Ishii et al [9] . 
Towards the end of section 2, we will lay down all the approximations that we 
have made, and also comment on some of existing important results in the area. 
Then, in section 3, we present our main numerical results. Section 4 ends this 
paper with discussions and directions for future research. 

2 Computation of Tidal Effects in Naked Sin¬ 
gularity Backgrounds 

We remind the reader that throughout this paper, we adopt geometrized units 
and set G = c = 1. Units will be explicitly restored when we discuss some 
examples towards the end of section 3. We will also use T(^ij)k to denote the part 
of the tensor T symmetric in indices i and j. Latin and Greek indices are used 
to denote spatial and space-time coordinates respectively. To begin with, we will 
illustrate the setup of the Fermi normal coordinates, following the work of |8]. 

2.1 Metric in Fermi Normal Coordinates 

In order to compute the tidal potential observed by a local inertial observer, we 
set up a locally flat system of coordinates on an arbitrary timelike geodesic, G. 
This would provide a way for a freely falling observer to report the effects of 
graviational held gradients in local experiments. We follow the standard proce¬ 
dure for constructing Fermi normal coordinates by setting up a tetrad basis at a 
point on the geodesic as the origin of our new coordinates, |pp where unprimed 
indices denote the Fermi normal coordinates. We parallely transport the tetrad 
basis along G and in particular choose eo(T) to be the tangent vector along G 
which is, by dehnition, parallel transported. We hence obtain the tetrad basis on 
the entire geodesic as a function of the proper time. The Fermi normal coordi¬ 
nates, x" of a point P' in the neighbourhood of the geodesic are then specihed. 
We dehne x° = r and x* corresponds to the point along the unique spacelike 
geodesic, G' at proper distance s with tangent vector at P{t) given by direction 
cosines x^/s in the tetrad basis, ep 

In order to demonstrate the differences in tidal effects, we will consider gen¬ 
eral static, spherically symmetric spacetimes to cover a broad class of physical 
scenarios. We present an explicit computation of the Fermi normal coordinates 
along circular and radial geodesics for these. The class of metrics that we consider 
in this paper are given by 

ds^ = —A{r)dt‘^ + B{r)dr‘^ + G(r)dr2^ , (2) 

where the coefficients are positive functions of the radial coordinate and dhl^ = 
d9^ + sin^ 9d(j)‘^ is the metric on the unit two sphere. 
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The metric possesses two cyclic coordinates, t and 0 corresponding to the 
conserved quantities 

E = A{r)i, L = C(r)(p , (3) 

where E and L are the energy and angular momentum per unit mass of the test 
particle. From the remaining geodesic equations, we reduce the problem to the 
equivalent one-dimensional problem with effective potential 


V{r) 


1 r ^ 

B{r) C(r) Ar 


(4) 


Imposing V{r) = V'{r) = 0 for circular orbits gives 


E = 




L = 




^A{r)C’{r)-C{r)A'{ry ^A{r)C'{r) - C{r)A'{r) ' 


(5) 


We now set up the Fermi normal basis which must satisfy the parallel transport 
conditions 

e^.e^ = j Cg' = 0 , (6) 

and Co is the tangent vector to the geodesic. This gives us the tetrad 


_ 

Cn — 


Cl = 


Co = 


63 = 


E L 

-,o,o. 


A{ry ’ ’^(r) 

Lsin(f2r) cos(r2r) 


VMr)Ciry ^/W) 

1 


, 0 , 


E sin(f2r) 
y/A{r)C{r) 


0 , 0 , 




:,0 


Lcos(f2r) sin(f2r) E cos (fir) 


^A{r)C{ry ’ ^A{r)C{r 

Here, the constant angular velocity of motion is 


Q = D{r)y = 


A'{r)C'{r) 
A{r)B{r)C{r) ’ 


where we have dehned 


D{r) = 


1 C'{r) [A{r)C'{r) - A'{r)C{r] 


A{r)B{r)C{r) 


In case of radial geodesics, proceeding similarly, we obtain 

E 


t = 


A{r) ’ 


r = 


'E^-A{r) 

A{r)B{r-) 


(7) 


( 8 ) 


(9) 


( 10 ) 
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We set up the tetrad satisfying Eq.(|6]) using the geodesic equations in this case 


eo' = 


Cl = 




= 0 , 0 , 




:,0 


= 0 , 0 , 0 , 


^/C{^ sin 9 


( 11 ) 


2.2 Calculating the Tidal Tensor 


Once the tetrad is constructed, we obtain the components of the curvature tensor 
in the tetrad basis 


Rafi'yS p'n' ■ (1^) 

The metric of an observer in the Fermi normal coordinate system can be expanded 
as 

1 1 ^ 


Qfiu T T q 9p,u,ijk^ ^ T 0{x ) . 


(13) 


The expressions for the derivatives of the metric in the Fermi normal basis were 
computed in [8] and [9]. The tidal potential, associated with the metric can 
then be expanded as 


= -^(S'oo + 1) = ]^Roiojx"x^ + ^RQ{i\o\j.k)x''x^x^ + 0{x'^) . (14) 


In case of the circular geodesic, the non-vanishing derivatives of the metric 
needed for the computation of the tidal tensor at leading order are 


S'oo.ii 

5'00,22 


CA'‘^C' - AC'^A' - 2ACC'A" + 2ACA'C" 
4ABC{AC' - A'C) 

A' ^ A'C 

A'C - AC' ^ ~ A ABC ' 


(15) 


For radial geodesic, the energy of the test particle F is a free parameter. This 
parameter cannot however assume large values, as we discuss in subsection 2.5. 
The corresponding non-zero quantities in this case are 


5'oo,ii 

^700,22 


A^[-2ABA" + A'{BA' + AB')] 

E‘^BCA'C' + A{E‘^ -1){BC'‘^ + C{B'C' -2BC")) 

900,33 — 4^2 Q2(J2 (^®) 


In our analysis, we will focus on the following class of static, spherically sym¬ 
metric metrics to compare the black hole and naked singularity cases- 
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• Janis Newmann Winicour (JNW) spacetime [13], [H]- This is a solution to 
the Einstein-Klein-Gordon system of equations. The spacetime consists of 
a singularity that is globally naked at r = B. 

ds^ = - (l--] dt^ + (l--] dr^ + r^(l--] dn^ . (17) 

\ r J V^/ V^/ 


This spacetime is sourced by the scalar field with magnitude q 






( 18 ) 


Here, 5 is a parameter related to the ADM mass by i? = and 

u = 2M/B lies between 0 and 1. Note that the case z/ = 1 corresponds to 
the Schwarzschild metric which is a black hole solution. For a fixed value 
of the ADM mass, moving away from v = 1 takes us deep into the naked 
singularity regime. Note however that there is a tradeoff here, namely that 
smaller values of v implies larger values of H, which in turn means that the 
solution is valid from a higher value of the radial coordinate r. Also, we 
would be mostly interested in stable circular orbits, and the criteria for such 
orbits are well established. We simply quote the result that defining the 
variable 

+ 3z/ ± \/5z/2 — Ij , (19) 

it can be checked that for 0 < z/ < l/\/5, stable circular orbits exits for all 
values of the radius (greater than B). For z/ lying between l/\/5 and 1/2, 
stable circular orbits are possible for radii less than r_ and greater than r+, 
while for v greater than 1/2, such orbits are possible for all values of the 
radius greater than r+. In our analysis on the JNW spacetimes, we have 
chosen values of the radii that correspond to stable circular orbits. 

• Bertrand spacetimes (BSTs) [T5|, [IS]- These are seeded by matter that can 
be given an effective two-fluid description and are a candidate for galactic 
dark matter. We look at the type II BST 

ds^ = - ^ ^ dt^ ^dr^ -F . (20) 

D + - r 


Here, is a rational. The parameters a and (3 can be used to provide a 
phenomenological estimate for the total dark matter mass of the galaxy (in 
good agreement of known results for low surface brightness galaxies) as [T7] 

aRl 

M = -^- , 

2(a + DRgf ’ 



where a good estimate for the size of the galaxy is Rg = a/D. It can be 
checked that here we have a central singularity at r = 0, which is naked. For 
computations performed with BSTs, we have set M = 1, which implies, via 
Eq. fl2l|) that a = 8D^. This will be understood in sequel. It can be checked 
that for BSTs, all circular orbits are stable. 
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• Joshi Malafarina Narayan (JMN) spacetime [18], [19] - describes a class of 
geometries that result from the gravitational collapse of a matter cloud with 
regular initial conditions into an astympotically static equilibrium conhgu- 
ration containing a central naked singularity at r = 0 

Mn 

ds^ = -(1 - ^o) ° ^ ^ 

This metric also matches smoothly to a Schwarzschild exterior at r = i?o 
with mass M = MqRq/2. The system possesses stable circular orbits for 
Mo < 2/3. 

While the JNW, BST and JMN space-times discussed above contain naked singu¬ 
larities, towards the end of this paper, we will also briefly comment on the nature 
of tidal forces in the interior Schwarzschild solution. This describes the metric 
of a fluid of constant density that is matched with an external Schwarzschild 
solution of mass M, and reads 


, 2 A 2Mr2\ \ 2 2 

df -F f 1- j dr^ + r^dQ . 

(23) 

The matching radius here is r = i?o- 


ds2 = 



2M 

Rn 



2Mr2 \ 


2.3 Hydrostatic equations for Computing Equilibrium Con¬ 
figurations 

We will now review a formulation to compute the equilibrium conhgurations 
under gravitational tidal forces arising from the potential gradients in a curved 
spacetime geometry. The discussion in this and the next subsection essentially 
follows [9] , and we closely follow the notation of that paper. We assume that the 
radius of the star El is much smaller than the radius of its orbit so that in Fermi 
normal coordinates, the self gravity of the star is described by Newtonian gravity, 
with an additional tidal potential due to the curved background. In such a case, 
the fluid star obeys the hydrodynamic equation 


dvi 


dxi 


dP 



( dAj 



dx^ 

1 

+ 


y dx^ 

dx^ / 

dr _ 


( 24 ) 


This equation is analogous to the Euler equation with the tidal potential 0* su¬ 
perposed with the Newtonian potential 0 and an additional term associated with 
the gravitomagnetic force described by the vector potential = ^Rk{ij)oX^xA 
Here, n* = dx^/dr is the fluid three-velocity and P is the fluid pressure. The 
Newtonian potential produced due to the mass of the star obeys the Poisson 
equation A0 = 47rp. 

^By “radius” of a stellar object, we really mean the length of its major axis. This slight 
abuse of notation should be kept in mind. 
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As is common in the literature, we assume a polytropic equation of state for 
the star given by 

P = , (25) 

where k and n are polytropic constants. The equation of continuity, along with 
that for hydrostatic equilibrium then leads to the Lane-Emden equation 


1 d f^2d0\ 


+ e^ 


0 , 


(26) 


where is a dimensionless radial parameter obtained by scaling the radius, and 
P = Pod"', where pc is the central density. In our computations, this will be used 
to £x the stellar radius in terms of the Lane-Emden coordinate at the stellar 
surface 


R 


[n + \)nfP ^ 

-te- ) 


(27) 


We will consider the values of the polytropic index n = 0.5, 1 and 1.5 for which 
it can be checked that .^o = 2.75, tt and 3.65 respectively. 

We hrst consider the case of circular geodesics. Following [9], we assume a 
co-rotational velocity held for the huid star with the velocity held 

u* = [—— a:csin(f2r), 0, — Xc cos(f2r)] , (28) 


where Xc is a small correction constant. In order to simplify numerical compu¬ 
tations, we go to a frame where the Euler equation is independent of r. This is 
achieved by choosing the frame, x® rotating at angular velocity about x^ axis. 
The hydrostatic equation corresponding to the hrst integral of the Euler equation 
in this frame is then given as 


— [(x^ — 2xc)^ -|- (a^^)^] — ^(ar -l- 1)^" + 0 + 0* + <(>m + C . (29) 


Here 0^ is the contribution due to the gravitomagnetic force and we have used 
the polytropic equation. 

In case of radial geodesics, El we take the velocity held to be irrotational, and 
hence the gradient of a velocity potential, 

^ • (30) 

The hrst integral for (12T)) is then given as 


dll) \ d^l) d^l) 1 


+ 0 + 


d^l) 


dx^ 


+ C 


(31) 


Here, '0 is determined by solving the equation of continuity. 


pA'd + 


^ dip dp 
c?x* dx^ 


0 . 


(32) 


These, along with the Poisson equation are the governing equations for the radial 
case. 


^There is an assumption of instantaneous equilibrium here. See the discussion in subsection 

2.5. 
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2.4 Numerical Procedure 

For our numerical procedure, we again closely follow |9]. We first consider the 
circular case. Here the governing equations are the Poisson equation and Eq. fl2^ . 
and these are to be solved iteratively. In order to achieve a convergence in the 
iteration, we switch to a dimensionless coordin_ate, x® = hg®. We correspond¬ 
ingly dehne the rescaled potentials, cf) = 0/h^, 0* = (l)t/h? so that our system of 
equations reduces to 

Ah(f> = 47rp , (33) 

02 _ _ _ 

— — 2gc)^ -|- (g^)^] = K{n + l)p" -|- + 4>t-\- <Pm) + C . (34) 

This system contains three primary constants, namely h, qc and C. The re¬ 
maining constants are hxed by choosing units so that the Newtonian mass of the 
central body is unity or by choosing the metric parameters in these units. The 
equilibrium conhgurations will be computed for different values of Pc- The three 
primary constants will have to be hxed at each step of the iteration. As in [9], 
we solve the Poisson equation in Cartesian coordinates on a uniform grid of size 
(2A^-|-1, A^-|-l, 2A^-|-1) in order to cover the region —L < g^ < T, 0 < g^ < L, 
—L < < L. Here, we have assumed rehection symmetry with respect to the 

g^ = 0 plane. We typically choose = 50 and the grid spacing to be gs/40, 
where g^ is the coordinate at the stellar surface along the g^ axis. We use a For¬ 
tran subroutine to solve the fourth-order hnite difference approximations to the 
elliptic partial difference equations. We implement a numerical algorithm that is 
summarized as : 

• We take a trial density distribution and use the cubic grid-based Poisson 
solver to solve the Poisson equation. We impose Neumann boundary condi¬ 
tions on the g^ = 0 plane to achieve rehection symmetry along the g^ axis 
and Dirichlet boundary conditions elsewhere, namely, 

0 —- [pd^r + 0{r~^) . (35) 

r J 

The distribution 0 is obtained which is used at the next step. 

• The constants h, g^, C are determined from the Euler equation by imposing 
constraints on the density prohle. We match the central density, p(0, 0, 0) = 
pc and require that dp/dq^ = 0 at the center. Also, we require the density 
to vanish at the stellar surface, p(gs,0,0) = 0. The set of constraints when 
substituted in (|3T|l can be solved to determine the values of the constants. 
From these values, the new distribution p(g*) is computed again using (|3T)l . 

• The new density distribution is truncated at the hrst instance where it goes 
to zero and is again used as a source in the Poisson equation. The process 
is repeated iteratively until it converges to the equilibrium distribution and 
the density prohle becomes stationary upto a numerical tolerance value. 

• We compute several equilibrium conhgurations for decreasing pc and track 
the quantity y; = dp/dq^ at the stellar surface (g^, 0, 0). The critical value of 
Pc corresponding to the tidal disruption limit i.e. the Roche limit is obtained 
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when X goes abruptly from a positive value to zero. Beyond this point, the 
density distribution becomes flat, signalling tidal disruption. For positive y, 
the star is in a stable conhguration. 

We can adopt a similar procedure for radial geodesics where there is an extra 
function that has to be determined. This can be handled using the additional 
constraint imposed by the continuity equation. We proceed as follows. As be¬ 
fore, a trial density function is used to solve the Poisson equation and obtain 0. 
This initial density distribution is also used to solve the elliptic partial differen¬ 
tial Eq. fl3^ using a three dimensional Cartesian PDF solver based on multigrid 
iteration. Once 0 and are known, we use Eq. (l3T]) to determine the contants h 
and C. This is done by imposing two constraints, namely p = pc at the center 
and the density must vanish at the stellar surface. One can now use the values of 
these cosnstants in Eq. fldTll to determine p on the entire grid. This new density 
distribution in again fed into the Poisson and continuity equations and the iter¬ 
ative procedure is used as previously to compute the equilibrium conhguration, 
and hence the Roche limit. 

2.5 Approximations, Limitations and Related Issues 

Before we move on to present our results, we point out the approximations that 
we will make, and the limitations of our analysis and possible caveats that need 
to be kept in mind. 

Since we are mainly interested in demonstrating the difference in tidal forces 
between naked singularities and black hole backgrounds, it suffices to work up 
to second order in the tidal potential. Fourth order effects can be included in 
our analysis, but these, apart from being small, will not qualitatively change our 
results. A more rehned analysis of the scenarios presented here including higher 
order corrections to the metric in Fermi normal coordinates is left for the fu¬ 
ture. Further, the contribution of the gravitomagnetic term in the hydrodynamic 
equations is only signihcant at fourth order in the expansion and so, these terms 
can be neglected at our level of approximation. Also, for the radial geodesics, for 
simplicity, we have chosen the velocity potential to be zero to give the irrotational 
held. In the rest of the paper, we will proceed with these approximations. 

There are a few caveats that we need to keep in mind. First, note that we 
work in a probe approximation, where the effects of the star back reacting on 
the metric is ignored. Effects of back-reaction are somewhat intractable in the 
present formalism and pose a formidable challenge. Here we proceed with the 
assumption that the background metric is hxed. 

Second, in case of radial geodesics, our analysis for the Roche radius is valid 
under the approximation of instantaneous equilibrium. In general, this radius will 
be less (or equivalently, the critical mass will be smaller) than that calculated 
from the numerical simulation by a small correction equal to the distance travelled 
by the star in the delay period during which the star attains the equilibrium 
conhguration. The approximation is valid in the limit that the distance travelled 
in the time scale during which the star attains equilibrium is smaller than the 
distance over which the tidal force changes appreciably. This condition is most 
naturally satished for low energy scales. We will, for radial geodesics, set E = 
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10“^ and it is assumed that at our level of approximation, the tidal force is 
constant over the range of validity of the Fermi normal expansion. 

Next, we note that a general tidal model can also result in a weak tidal 
encounter between a star orbiting in a massive background. Such encounters can 
result in a mass loss from the star, as discussed in details in [20]. These effects on 
white dwarfs are seen to increase the likelihood of tidal disruption as a function of 
time. In this work, we will however restrict ourselves to conditions for complete 
tidal stripping for both radial and circular orbital motion. 

Finally, although not presented here, the present analysis can be used to build 
upon the work of izu to calculate rates for tidal disruption events for supermas- 
sive black holes capturing not only the relativistic treatment of tidal disruption 
but also the hydrodynamics of the orbiting star. The method is summarized as 
follows. Monte Carlo simulations are performed with an appropriate distribution 
of the free variables in the geodesic equations for a general orbit. These equa¬ 
tions, in Fermi normal coordinates, are integrated upto the pericenter. In our 
case, we can run the numerical code to determine if the star is directly captured 
or is tidally disrupted at this point by computing the equilibrium hydrodynamic 
conhguration. This should give us an estimate for the tidal disruption rate of 
the simulated orbits. One can then study the effects of varying different metric 
parameters and astronomical data from flares associated with these tidal disrup¬ 
tion events can be used to provide insight into the nature of the background 
singularity. In principle, this should improve upon the work of ED. but such an 
analysis is left for the future. 


3 Results and Analysis 

In this section, we will present the results of our numerical analysis. We remind 
the reader that unless mentioned otherwise, in all examples below, we have set 
the radius of the star to be 0.5, and the mass of the background to unity. We 
start with circular geodesics in various space-times. In Fig.([2|), we show the 
critical mass as a function of the radial distance for JNW space-times, with the 
polytropic index (of Eq. (l25|l ) n = 0.5. The solid red, dotted green and dot-dashed 
blue curves correspond to the values z/ = 0.01, 0.1 and 0.6 respectively, in Eq. (lT7|l . 
For comparison, we have shown with the dashed black curve the corresponding 
result for the Schwarzschild black hole (z/ = 1 in Eq. (11711 h Clearly, the effect 
of this naked singularity background is seen to increase the tidal force on the 
star. This is our Erst observation : if, at a given radius, stellar objects above the 
critical mass predicted from the Schwarzschild black hole are seen to exist, they 
might point towards a naked singularity, rather than one with a horizon. 

In Fig. (El), the computation is repeated for the BST naked singularity back¬ 
ground of Eq. fl20|) . In this hgure, the solid red, the dotted green and the dot- 
dashed blue curves correspond to choosing the values (3 of Eq. fl20|) as 0.2, 0.5 
and 0.8, respectively. We see here that increase in the value of (3 corresponds 
to higher tidal forces. This is again indicative of the fact that for similar cen¬ 
tral masses, the BST naked singularity predicts higher tidal disruption limits. If 
observational indications of this fact are found in future experiments, BSTs can 
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Figure 2: Critical Mass for circular 
geodesics in JNW space-times for differ¬ 
ent u, with n = 0.5 



Figure 4: Leading Order Tidal Potential 
for Circular BST, with the chosen values 
D = 10^ 13 = 0.8, r = 3. 


Figure 3: Critical Mass for circular 
geodesics in BSTs for different f3, with 
D = 10^ n = 0.5 



Figure 5: Leading Order Tidal Potential 
for Circular BST, with the chosen values 
D = 10^ 13 = 0.8, r = 8. 


possibly be used as models to understand such a scenario. 

Next, we come to the tidal potentials, computed from Eq. ffTdll . We show this 
for two examples. Figs.(jl]) and (|5]) shows the variation of the second order tidal 
potential along the three axes in the Fermi normal coordinates, for r = 3 and 
r = 8. Although the plots have been shown for BST, we have checked that the 
nature of the tidal potential is similar for other metrics, for circular geodesics. 
The gradient of the plot is a measure of the tidal force. The given plots indicate 
that the potential is conhning along the and axes and hence the resulting 
equilibrium conhguration will have major axis along q^. Note that the leading 
order tidal potential has a symmetry about the origin along all the axes. This 
symmetry is generically broken in small amounts when higher order terms in 
the tidal potential are taken into account. We observe from Figures ([1]) and ([5]) 
that as r is decreased, the tidal potential becomes steeper, indicating higher tidal 
effects. For the Kerr case, we observed that the potential is less steep along the 
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Figure 6: Density variation near equi- Figure 7: Density variation near equi¬ 
librium for BST (Circular Orbit), with librium for BST (Circular Orbit), with 
D = 10^ (3 = 0.8, n = 0.5 D = 10^ [3 = 0.8, n = 1.5 



Figure 8: Critical mass for radial Figure 9: Critical mass for the radial 
geodesics in the JNW background for geodesics in the JMN background for dif- 
different polytropic index, where we have ferent values of the radius of the star Rq 
taken z/ = 0.4. with Mq = 0.1 

and axes for a higher value of a indicating that that the tidal forces weaken 
as the effect of spin is taken into account. It is also interesting to look at the 
behavior of the density variation of the star near equilibrium. Figs.(l6]) and fl71) 
show the typical variation of the equilibrium density distribution at the Roche 
limit in the second order tidal potential along the q^ axis for different values of 
where we have chosen the polytropic index to be 0.5 and 1.5, respectively, as 
a function of the coordinate q^. In Fig. ([6]), this is shown for = 0 (solid 

red), g^ = 0,g^ = 3R/8 (dashed green), g^ = 0,g^ = R/2 (dot-dashed blue) and 
g^ = 0,g^ = 5i?/8 (sparse-dotted magenta). In Fig.([7]), the same values of g^ 
and g^ are used for the solid red and the dashed green curve. In this hgure, the 
dot-dashed blue curve corresponds to g^ = 0, g^ = R/4. 

The density in equilibrium is set to zero at g^ = 0.8 along the major axis on 
the grid corresponding to Rq. As expected, the density along the g^ direction goes 
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Table 1: Comparative Table For Circular Geodesics for r = 8 


Metric 

Parameters 

Merit 

Schwarzschild 

M = 1 

2.5 X 10-3 

JNW 

o 

bo 

2.8 X 10-3 

JNW 

V = 0.6 

3.5 X 10-3 

JNW 

1/ = 0.4 

6.1 X 10-3 

BST 

D = 10^, 13 = 0.1 

5.7 X 10-^ 

BST 

D = 10^, 13 = 0.6 

1.9 X 10-3 

BST 

D = 103, (3 = 0.8 

5.4 X 10-3 

BST 

D = 10^ (3 = 0.3 

2.8 X 10-3 

JMN 

Mo = 0.2 

1.1 X 10-3 

JMN 

Mo = 0.08 

4.4 X 10-^ 

JMN 

Mo = 0.02 

1.1 X 10-^ 


to zero faster than in the direction. We observe that the density is symmetric 
with respect to the = 0 plane. This symmetry is generically broken slightly 
when higher order terms in the tidal potential are taken into consideration. For 
smaller values of the polytropic index, the pressure, in accordance with the stiffer 
polytropic equation of state, has a stronger density dependence. This is reflected 
in the fact the equilibrium conhguration is more compact with higher surface 
density for n = 0.5 whereas for n = 1.5, the density goes to zero smoothly at the 
surface due to lower pressures. 

Finally, in Figs.(|H]) and ([9]), we show the critical mass versus radial distance for 
the JNW and JMN singularities. In Fig.®, we show the variation of the critical 
mass for different values of the polytropic index, with the solid red, dotted green 
and dot-dashed blue curves corresponding to n = 0.5, 1 and 1.5 respectively. In 
Fig.®, this is shown as a function of the radius of the star, with the solid red, 
dashed green and dot-dashed blue lines corresponding to this radius being 0.5, 
0.75 and 1, respectively. 

We now tabulate some numerical data on the various cases that we have dis¬ 
cussed. In table CDi, we have provided such data for circular geodesics in the 
Schwarzschild, JNW, BST and JMN metrics, with the mass of the singularity 
taken to be unity in all the cases. One can see from the table that at this value 
of the radial distance, the tidal force can be of similar orders of magnitude for 
the Schwarzschild black hole, and the JNW and the BST naked singularity back¬ 
grounds. Specihcally, tidal forces are seen to be stronger as one decreases the 
value of V in the JNW background. For example, aX v = 0.4, the critical mass al¬ 
lowed by this background is about three times that in a same-mass Schwarzschild 
background. This would imply that if a stellar object of higher critical mass than 
predicted by a Schwarzschild analysis is found, this might indicate that one needs 
to look closer at the nature of the central singularity. On the other hand, the JMN 
background indicates that stellar objects having masses much smaller (by even 
two orders of magnitude) than those predicted by a Schwarzschild analysis might 
exist in stable circular orbits. Although direct observational evidence for such 
objects might be extremely difficult in practise, nonetheless our simple minded 
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Table 2: Comparative Table For Radial Geodesics for r = 8, E = 0.01 


Metric 

Parameters 

Merit 

Schwarzschild 

M = 1 

1.5 X 10-3 

JNW 

o 

bo 

1.9 X 10-3 

JNW 

v = 0.6 

2.1 X 10-3 

JNW 

1/ = 0.4 

4.7 X 10-3 

BST 

D = 10^, 13 = 0.1 

1.4 X 10-^ 

BST 

D = 10^, 13 = 0.6 

3.4 X 10-3 

BST 

D = 103, (3 = 0.8 

3.9 X 10-3 

BST 

D = 10^ (3 = 0.3 

1.3 X 10-2 

JMN 

Mo = 0.2 

5.4 X 10-^ 

JMN 

Mo = 0.08 

2.7 X 10-^ 

JMN 

Mo = 0.02 

8.1 X 10-3 


analysis indicates nsefnl bonnds on the masses of stellar objects in circnlar orbits 
of a given radins. Of conrse, we have conhned onr analysis to circnlar orbits, 
whereas actnal orbits of stellar objects might be highly elliptical. We however 
expect that qnalitative featnres of onr analysis will remain nnchanged in snch 
sitnations also, thongh this reqnires a more detailed analysis. 

As a conseqnence of onr analysis, we note that for large valnes of r, the 
differences in the magnitnde of the tidal force between the Schwarzschild and 
the JNW naked singnlarity backgronnd increases appreciably. For example, for 
circnlar geodesics at r = 12, we hnd that whereas Merit = 9.8 x 10“^ for the 
Schwarzschild case, it is 5.6 x 10“^ for the JNW naked singnlarity with u = 0.2. 
For circnlar geodesics at r = 21, the Schwarzschild critical mass is 8.6 x 10“^, while 
the JNW backgronnd in this case with i/ = 0.1 predicts a critical mass 5.5 x 10“^, 
almost an order of magnitnde difference. This latter set of nnmbers translate 
into interesting realistic ones as we now illnstrate. Take the mass of the central 
singularity as 1.8 x IO^Mq, and Rq = 14 km (which correspond to i?o = 5 x 10“^ 
in natural units). Then, r = 21 in natural units leads to r = 5.6 x 10^ km. The 
Schwarzschild background in this case predicts a critical mass Merit = I.OMq, 
while a JNW background with i/ = 0.1 predicts Mcnt = 15.7Mq, where we have 
assumed n = 1 in the polytropic equation of state of the neutron star. Given 
a typical neutron star mass ~ 1.5M0, we see that the neutron star is tidally 
disrupted for the naked singularity background but not for the Schwarzschild 
background at this radius. 

It is also worthwhile to mention that while the discussion of the previous 
paragraph focused on the JNW naked singularity, the result assumes signihcance 
given the fact that similar conclusions can be reached for BSTs as well. Indeed, 
from Fig. ([3]), we see indications that the critical mass increases sharply as the 
value of the parameter (3 in Eq. (l20|) is increased. It is thus expected that appro¬ 
priately tuning the parameters of a BST, one can in principle obtain values of the 
critical mass substantially larger than what is obtained in Schwarzschild back¬ 
grounds. On the other hand, as Table (P) indicates, JMN backgrounds typically 
show lower tidal effects compared to the Schwarzschild cases. 
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Table ([2]) summarizes the results for the radial geodesics, which have to be 
understood along with the limitations discussed in subsection 2.5. We see that 
these are qualitatively similar to those for the circular geodesics. In this table, 
we have considered a small energy of the star. Increase in this leads to a higher 
difference in magnitude for the critical mass for the naked singularity background 
compared to the Schwarzschild case. However, the assumption that the star has 
travelled a distance that is small compared to the scale over which the tidal forces 
change will not be valid for high energies of the star, and thus higher values of 
energy will probably not be very trustable in our framework. 

Before we end, let us briefly comment on the nature of tidal forces in the 
interior Schwarzschild solution of Eq.(|23]). We considered radial geodesics in this 
geometry, with energy E = 0.01. Here, we have taken a matching radius r = 20. 
Our results in this case suggest that the critical mass reaches a saturation value 
as the radial distance of the stellar object decreases. However, we were able to 
achieve good numerical convergence in this case only for a limited number of 
parameter values. This case is of interest as this exemplihes motion in a diffused 
matter background by a stellar object (which is assumed not to back react on the 
background) and merits further understanding, which we defer for a later study. 


4 Discussions and Conclusions 

In this paper, we have examined the nature of tidal forces in a class of non-rotating 
naked singularity backgrounds on stellar objects in radial and circular geodesic 
motion. The purpose of this work was to understand theoretical differences in the 
nature of these forces, keeping in mind observational aspects. Broadly, this paper 
establishes the magnitudes of these forces for three different naked singularity 
backgrounds. Our main conclusion here is that tidal forces can be signihcantly 
different for these cases when compared to the Schwarzschild background. If in 
experiments can indicate numerical bounds on the masses of objects that can be 
tidally disrupted by a singularity, our results might be effective in modelling such 
situations, to indicate the nature of the underlying singularity. 

In this paper, we have conhned ourselves to a second order expansion of the 
tidal potential in Fermi normal coordinates. The analysis of [9] , which is based on 
a fourth order expansion of the potential is numerically more accurate. However, 
in this paper our main focus was comparing the magnitudes of the tidal force in 
different backgrounds. Inclusion of higher order corrections, although important, 
will not qualitatively affect the results of the present analysis. 

Also, our analysis is simple minded, has a number of constraints which we 
have extensively discussed. Nevertheless, we believe that the present analysis 
complements the work of [9] and sets the stage for a deeper question, namely if 
tidal disruptions can be an effective indicator of the central singularity in galaxies 
or galaxy clusters. As an immediate future application, we hope to build upon 
the results of this paper to understand rates of tidal disruption events, while 
taking into account the fluid nature of the star. This work is in progress. Fur¬ 
ther, as is well known, a realistic neutron star has a superfluid character, and its 
fluid dynamics might be different from what is assumed in this work. How this 
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will affect tidal forces on neutron stars might be an important issue to understand. 

Note added : A preliminary version of the FORTRAN code used for the com¬ 
putations in this paper is available upon request. 
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